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ABSTRACT 

The plasmoid-induced-reconnection model explaining solar flares based on 
bursty reconnection produced by an ejecting plasmoid suggests a possible relation 
between the ejection velocity of a plasmoid and the rate of magnetic reconnec- 
tion. In this study, we focus on the quantitative description of this relation. We 
performed magnetohydrodynamic (MHD) simulations of solar flares by changing 
the values of resistivity and the plasmoid velocity. The plasmoid velocity has 
been changed by applying an additional force to the plasmoid to see how the 
plasmoid velocity affects the reconnection rate. An important result is that the 
reconnection rate has a positive correlation with the plasmoid velocity, which 
is consistent with the plasmoid-induced-reconnection model for solar flares. We 
also discuss an observational result supporting this positive correlation. 

Subject headings: MHD — Sun: corona — Sun: flares 



1. Introduction 



Magnetic reconnection is a process in which the magnetic energy is converted into the 
kinetic and thermal energy (ISweetl Il958l ; iParkerl Il963l ; iPetschekl Il964l ). and it has been 
widely believed to play a fundamental role in causing solar flares. A model for flares 
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based on magn e tic reconn e ction has b een developed sinc e 1960 s by ICarmichaell (119641 ). 
Sturrockl (119661 ). iHirayamal (119741 ). and iKopp fc Pneumanl (119761 ). so this model has been 



called the CS HKP model. The ob servations supporting t 



such as cusps (ITsuneta et al.lll992al ). arcades (ITsuneta et al. 



lis model has been reported , 



1992 



b; 



Isobe et al. I2OO2I ), loop top hard X-ray (HXR) sources ( Masuda et al 



McAllister et al 



19941 ; 



1992; 



Sui fc Holman 



20031). X-ray jets (IShibata et al.lll992l ; IShimojo et al.lll996l ). and so on. Furthermore Yohkoh 
( lOgawara et al.lll99ll ) discovered the common property of flares with diff erent appearan ces, 
such as Long Du ration Events (LDE flares) and impulsive flares (see IShibatal Il999l and 
Aschwanden 2002 for review). 



So far the dynamic process ca used by magnetic reconnection in fl a res has been widely 



studied by MHD simulations 


e.g. 


Forbes & Priest 


1983: 


Forbes 


1990l: 


Forbes & Malherbe 


1991; 


Maeara et al. 


1996; 


Ueai 


199d; 


Lin & Forbes! 


2000; 


Yokovama & Shibata 


2OO1I). Most 



of these works have been focused on MHD processes producing apparent features of flares 
such as cusp shaped loops, ejecting plasmoid (blob of plasma), inflows and loop-top HXR 
sources. However, these works have not clearly explained the fundamental physical process, 
that is what determines the energy release rate in flares. This question is related to the rate of 
magnetic reconnectio n, so to identify the condition under which fast magnetic reconnection 
(jUgai fc Tsudalll9771 ) operates is important. 



Observationally, solar flares are often associated with plasmo id ejections. IShibata et al. 



(119951 ) found that 8 impulsive flares on the limb (Masud a -type: iMasuda et al 



19951 ) were 



associated with plasmoid ejections. lOhyama fc Shibatal (119971 . Il998l ) carefully analyzed 
plasmoid ejections in impulsive flares and found that plasmoi d s und ergo strong accelera- 
tion during the impulsive phase of these flares. ISakajiri et al.l (120041 ) found the tiny two- 
ribbon flare driven by emerging flux accompanying the m iniature filament eruption ^plas- 



moid ) . This feature was also f ound in other observations flZhang et al 



2001 



Takasaki et al. 

20041 ; ISterling fc Moord 120051 ) and numerical simulations ( Magara et al.l Il997t ) . It is also 
found that t here is a positive correlation between the p l asmoid velocity and the recon- 



nection rate (IShibata et al.l Il995l ; iQiu fc Yurchvshvn 



other literatures discussing this topic 



Priest fc Forbes 2002; Lin et al. 2003 



2005 



Nitta &; Akivama 



Kim et al. 



20051 



199 



Forbes 



2000; 



Lin Sz Forbes! (120001 ) derived an an- 



Shimizu et al.l 120081) . There are 

IT 



Klimchuk 2001 



alytic relation between the acceleration of coronal mass ejections (CMEs) and reconnection 
rate. 

The observations above show an important suggestion that plasmoid ejection plays 
a key ro l e in c ausin g fast magnetic reconnection. Based on these observational results, 
Shibatal (119961 . Il997l ) extended the classical CSHKP model and proposed the plasmoid- 
induced reconnection model. In this model, a plasmoid (or flux rope in a 3D situation) 
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is created in the anti-parallel magnetic field by the magnetic reconnection (Figure [T^,, b, c). 
Then the plasmoid situates in a current sheet inhibits inflows into the sheet, so reconnection 
is inefficient and magnetic energy is stored (Figure dti). Then the plasmoid starts to move 
at the velocity v p i asmoi( i, inflows toward the X-point (v in fi ow ) are induced following mass 
conservation and reconnection starts (Figure [Tfe). If we assume the incompressibility, the 
mass flux into the reconnection region is given by ~ Vi n fi ow L in fi ow and the mass flux ejected 
by the plasmoid is ~ v plasmoid W p i asmoid , and these are balanced, where W plasmoid is the typical 
width of the plasmoid, and L in fi ow (> W p i asmoi d) is the typical length of the inflow region. 
Consequently the induced inflow speed can be estimated as follows: 



"^inflow ~ Vplasmoidft^plasmoid./ Linflc 



T 



Since the reconnection rate is determined by the speed of the inflows, fast reconnection 
becomes possible when plasmoid ejection is fast. Moreover the jet from a reconnection point 
accelerate the plasmoid □, so the fast reconnection further drives fast plasmoid ejection. This 
suggests a positive correlation between the plasmoid velocity and the reconnection rate. The 
merit of this model is to provide us with a unified view for understanding various types of 
flares with different spa tial sizes and timescales, and it can naturally cause fast reconnection 
(IPriest fc Forbes! Il986l ) in a c urrent sheet where many plasmoids with different sizes a re 
created by tearing instability (jKliem et al.ll2000l ; IShibata fc Tanumall200ll ; iKarlickyl 120041 ) . 



So far there have been no MHD simulations performed to examine the plasmoid-induced- 
reconnection model, which should reproduce the correlation between the plasmoid velocity 
and the reconnection rate. In this study, we performed a series of these MHD simulations 
by changing the parameters related to resistivity and plasmoid velocity to investigate the 
relation between these two quantities. 

The model and the numerical method are described in section [21 and the numerical 
results are presented in section [3j The discussion is given in section [4] and the conclusions 
are given in section [51 



^^Note that in addition to the acceleration by the reconnection jet, the magnetic pressure gradient force 
can also accelerate the plasmoid if there is a global magnetic pressure gradient around the plasmoid. See 
more detailed discussion in Appendix. 
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Plasmoid 



Fig. 1. — Schematic picture of the plasmoid- induced reconnection model. Solid lines indicate 
magnetic field lines. Panel (a), (b) and (c) shows the process creating the plasmoid in the 
anti-parallel magnetic field by the magnetic reconnection. Panel (d) shows how the plasmoid 
in the current sheet inhibits the magnetic reconnection in the reconnection region shown by a 
light gray area in panel (e). Thick black arrows indicate inflows to a reconnection point, and 
thick white arrows indicate reconnection jets. Panel (e) explains how the plasmoid ejection 
can induce strong inflows into the reconnection region. 
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Fig. 1. — cont. 
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2. Numerical Method and Model 
2.1. Numerical Method 



We performed simulations with a multistep implicit scheme (jHulll989l ) in order to solve 
the 2. 5- dimensional time-dependent resistive MHD equations in the Cartesian coordinates 
(x, y, z) where y is directed upward: 

|| + V.(pv) = 0, (2) 
<9v 1 1 

— + (v ■ V)v + -Vp j x B — F = 0, (3) 

at p p 

<9B 

— + V x (v x B) - 77V x j = 0, (4) 
at 

^+v-VT+ 7 -lTV-v-^ %j = 0, (5) 

at p 

where p is the gas pressure, B = V x (ipe z ) + (0, 0, B z ) is the magnetic field (ip is the 
magnetic flux function), j = V x B is the current density, rj is the resistivity, and F is an 
additional force which is introduced in section 12.41 The seven independent variables are the 
plasma density (p), temperature (T), velocity (v x , v y , v z ), magnetic flux function (ip), and 
perpendicular component of magnetic field (B z ). 

The units of length, density and temperature are L = 1 x 10 9 cm, p = 1.6 x 10~ 15 g 
cm -3 (i.e., Uq = 1 x 10 9 cm -3 ), and T = 1 x 10 6 K, respectively. Velocity is normalized by 



v = C s = J^Ml = 128.6 km s" 1 (6) 
V m H 

where ks is Boltzmann constant and rriu is the mass of a hydrogen atom. 

The parameter /3q is the typical ratio of gas to magnetic pressures in the lower corona 
above an active region. j3 is chosen to be 0.01 measured at (x,y) = (0,0.12) (magnetic 
configuration is described in section l2~2l) . and then the unit of magnetic field strength is 



Bo= / 16 £ AMi = 25 78 G 
m H /3 
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The unit of time is given by tao = L /vao, where vao is the Alfven velocity (vao = 
B /s/4:7rpo = 1818 km s" 1 ) and tao = 5.5 s. The resistivity (77) is described in section 

The size of the simulation box is —8 < x < 8 in the horizontal direction and < y < 15 
in the vertical direction. The domain is divided by 201 x 601 grid points, which are uniformly 
distributed in the y-direction while they are not distributed uniformly in the x-direction. 
The grid spacing Ax ranges from 6.46 x 10~ 3 to 3.14 x 10" 1 with the finest interval around 
x = 0. At the bottom of the simulation box is applied a line-tying boundary condition 
where all quantities except T are fixed; T is determined by extrapolation. At the other three 
boundaries is applied an open boundary condition. 



2.2. Initial Condition 



We assumed an octapolar magn etic configuration as th e init ial condition (Figu re [2]). 
Our model is based on the model of IChen fc Shibatal fl2000i ) and IShiota et all (120051 ) that 
assumes a quadrupolar configuration. We added two dipole fields to this model, in order to 
control the magnetic flux in the inflow region, duration and released energy of a flare. Short 
durations are the feature of impulsive flares. 

To get the magnetic configuration with a detached plasmoid (or a flux rope in 3D 
situation) and two side loops, our initial state consists of three separate groups of currents in 
a two-dimensional Cartesian plane: line current placed at the point (0, h) with finite radius 
r , its image current below the photosphere, and a background octapolar potential field 
produced by eight line currents below the photosphere, i.e., at eight positions (±0.1, —0.3), 
(±0.6,-0.3), (±2.0,-0.3), and (±2.5,-0.3), with the direction being "+, positive z", "-, 
negative z" , "+", "±", and "+" in sequence. Then, the initial magnetic 

configuration is expressed as 

rl> = ip b + il>i + ipi, (8) 

where the background field (ipb) and the magnetic components of the image current (ipi) and 
of the line current ("0z) have the following forms: 



06 



Ci In 



[(x + 0.1) 2 + (y + 0.3) 2 ][(x - 0.1) 2 + (y + 0.3)' 



+ c 2 ln 



[(x + 0.6) 2 
[(x + 2.5) 2 



(y + 0.3) 2 ][(x 
(y + 0.3) 2 ][(x 



0.6) 2 + (y + 0.3) 2 ] 
2.5) 2 + (y + 0.3) 2 ] 



[(x + 2.0) 2 + (y + 0.3) 2 ][(x - 2.0) 2 + (y + 0.3) 2 ] 



-c 3 y ln[x 2 



(y + hy 



(9) 
(10) 




-4 



-2 





x 



Plasma beta: 



0.1 1 


10 











Fig. 2. — A portion of the initial magnetic field configuration is presented. The size of the 
simulation box is —8 < x < 8 in the horizontal direction and < y < 15 in the vertical 
direction. Solid lines show magnetic field lines. ® and denote the direction of currents. 
Gray scale shows the plasma beta, which is the ratio of gas pressure to magnetic pressure. 
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r 2 



ifc = { C4 W (r " ro) ' (ii) 

C4 [r /2 - r ln(r ) + r ln(r)] , (r > r ), 



where r = ^/x 2 + (y — h) 2 is the distance from the center of the plasmoid, h is the height 
of the plasmoid and r is the radius of the plasmoid, which determine the configuration of 
magnetic field. In this study, we set h = 1.0 and r = 0.2. 

The coefficients ci, c 2 , c 3 and C4 in the equations (jUJ), (IIP)) and ffTTj) control the strength 
of the magnetic field. When the radius of the plasmoid, r , is small en ough compared to the 



height of the plasmoid, h, equilibrium solutions are analytically given by lLin &: van Ballegooijen 



( 120051 ). However, it is difficult to determine the coefficients analytically in our model, be- 
cause the plasmoid is in the complex background magnetic field and r is not too small. 
Therefore we determine the coefficients, c\, C2, C3, and C4, by trial and error in order to make 
the plasmoid center keep roughly stable and move much slowly compared to a rapid motion 
driven by reconnection. In this study, we set the coefficients as c\ = 2, c 2 = 2, c 3 = 6 and 
c 4 = 2.4. 

The physical meaning is as follows. If we remove the plasmoid, the magnetic field 
configuration are potential octapolar field, produced by the eight line currents, "0&, and the 
image current below the photosphere, ipi. The configuration makes a null point, which is 
saddle of potential above the photosphere. If we place a plasmoid with zero radius at the 
null point, then the plasmoid does not begin to move itself without any perturbation. 

Although it is difficult to create a state of equilibrium when a plasmoid has a finite 
radius, this is not a major problem because the purpose of this study is not to investigate 
the trigger mechanism of a flare, but to focus on the physical processes working after magnetic 
reconnection starts. 

To satisfy the force balance within the plasmoid, a perpendicular magnetic component 
(i.e., B z ) is introduced inside the plasmoid: 



b, = < a2 V 2 ( 1_ S)' r " ro; < 12 > 

0, r > r , 

Other quantities are set to be uniform; (p, T, v x , v y , v z ) = (1.0, 1.0, 0.0, 0.0, 0.0). 
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2.3. Resistivity 

The resistivity, 77, is assumed as an anomalous resistivity: 

v d > 2v c , 

V = { Vo I ^ - 1J , v c <v d < 2v c , 




(13) 



where v d = \j z \/p is the (relative ion-electron) drift velocity, and v c is the critical velocity. 
Here we assumed v c = 0.5. The anomalous resistivity begins to work when the drift velocity 
exceeds the critical velocity, otherwise there is no resistivity, only small numerical resistivity 
remains. It is known tha t the anomalous resistivity may be caused by plas ma instabilities 
(ICoppi &: Friedlandlll97ll ). and localized resistivity triggers fast reconnection (jUgai fc Tsuda 
19771 ). We assumed that r/ is a free parameter (see Tabled]). We simulated with various 
r)o in order to change the reconnection rate (case A). Therefore we can examine how the 
reconnection rate influences the plasmoid velocity. We also set the resistivity to be zero in 
order to prevent magnetic reconnection and investigated the role of reconnection in plasmoid 
ejection (case C). 



2.4. Additional Force 



Considering the plasmoid-induced- reconnection model (jShibatal Il996l . 119971 ). we can 
expect that there is a positive correlation between the plasmoid velocity and reconnection 
rate. To see this, we changed the vertical velocity of the plasmoid by using the vertical 
additional force that acts only inside the plasmoid (or a flux rope in a 3D situation) in case 
B. There is no additional force in case A and case C. The additional force is introduced 
as a virtual force in order to change the plasmoid velocity and see how it relates to the 

Table 1. Free parameter 770 which determines the resistivity and strength of the additional 

force (F y ). 



Case 


Fy 


Vo 


A 





0.2, 0.02, 0.002, 0.0002 


B 


±0.5, ±1.0 


0.002 


C 
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reconnection rate. In the present study we do not focus on the role of the additional 
force in triggering a flare, rather we investigate how the force affects the evolution following 
the onset of a flare. In fact, even if there is no additional force applied to the initial state, 
the plasmoid begins to move upward due to the non-equilibrium of the initial state. The 
additional force F in equation Q is defined as follows: 

f (0,F y F ,0), r<r ; 

1 (0,0,0), r >r , 1 } 

where Fq = 0.271poLo/C s , and F y is given in Table [TJ Note that the value of F y is also 
set to be negative, which means that the direction of additional force is downward and it 
decelerates the plasmoid. 



3. Results 



In this section, we introduce the typical result from case A. Figures [3] and [4] show 
the magnetic field distribution on x-y plane and Figure [5] shows the temporal variation of 
physical values in case A with rj = 0.002, which is called the standard case here after. 
Because of the non-equilibrium of the initial state, the plasmoid begins to move upward a 
little in < t < 5 in Figure [5b, UJ, [8b, [9b (i.e. melon-seed effect). After that the plasmoid 
starts to be accelerated when the reconnection starts until t ~ 15. The rising plasmoid 
causes the inflow of magnetized plasma into the magnetic null point (X-point). Magnetic 
field lines at both sides of the null point are carried into the point, and a current sheet is 
formed near the X-point. The drift velocity inside the current sheet then increases quickly 
by the inflow in 2 < t < 5 (Figure [5b,). When the condition of equation (fl3j) is satisfied, 
anomalous resistivity sets in and fast reconnection occurs (Figure [5b). This accompanies 
Y-shape and inverted Y-shape slow shocks at both ends o f the current s heet (Figure [6]), 
which is one of the features of Petschek-type reconnection (jPetschekl Il964j ). The plasmoid 
i s accelerated rapidly in 5 < t < 15 (Figure Eb). This result is consistent with observations 



(lOhvama fc Shibatal Il997|, 



1998 



Takasaki et all 12004 ; IZhang et all l200ll ; ISterling fc Moore 



20051 ) and a simulation (IMagara et al.lll997l ). However the plasmoid is decelerated gradually 
during 15 < t < 40, although reconnection jets push the bottom of the plasmoid (the white 
line in Figure [7b). This deceleration is mainly due to gas pressure, since the plasma beta 



2 The additional force used in the present simulations is a purely virtual force, and we do not intend to 
model the force realistically. For the interested reader, however, it may be useful to note that there have 
bee n proposed seve ral ways of modeling the force in a 3D space, such as curvature force (see lShafranovlll966 



and lLin et al.lll998l ) 
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(b) t=10 



(c) t=15 




-1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0 
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Pressure: 



0.0 0.5 1.0 1.5 2.0 



Velocity: 



Fig. 3. — Pressure distribution in a typical case (case A with 7/0 = 0.002). The solid lines 
denote magnetic field lines, and arrows represent velocity field. 
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above the plasmoid is high (see plasma beta distribution in Figure [2]). The rising plasmoid 
induces a strong inflow into the X-point in 5 < t < 40 (Figure [5U), and consequently the 
reconnection rate grows up. Here we define the reconnection rate as flux canceling speed i.e. 
dip/dt(= rjj) at the X-point. 

Figure [7] shows the time evolution of pressure and vertical velocity at x = 0. The bipolar 
jets from X-point are created by reconnection. The velocity of the upward jet reaches up 
to the local Alfven velocity. The jets created a fast-mode shock when colliding with the 
magnetic loop and the plasmoid. The upper fast shock is weaker than the lower one, because 
the relative speed between the reconnection jet and preexisted plasma is smaller. The lower 
fast-mo de shock (terminatio n fast shock) may be observed as a loop top hard X-ray (HXR) 



source (IMasuda et al.l Il994l ). The termination fast shock is rising as the magnetic loop 



develops due to the piling-up of reconnected field. 

Figure [8] shows time variation of the reconnection rate and the plasmoid velocity in case 
A with different values of rjo- The reconnection rate and the plasmoid velocity are enhanced 
as t)q grows. The reconnection rate seems to be proportional to log(r/) except rj = 0.2, but 
a careful analysis with wider parameter survey is needed for a final conclusion. 

The time variations of the reconnection rate, the energy release rate, and the plasmoid 
velocity are shown in Figure El when we changed the plasmoid velocity due to the additional 
force (case B). An interesting point suggested by this figure is that the plasmoid velocity 
changes significantly when the additional force is applied, but the reconnection rate does 
not. On the other hand, the energy release rate (Bf n j- low /4:7i x Vi n fi ow L in fi ow near the current 
sheet) is changed significantly, because the length of the current sheet (L in fi ow ) is changed 
by the position of the plasmoid. 

Figure [Sb and Figure EH show the correlation between the reconnection rate and the 
plasmoid velocity measured when the reconnection rate becomes the maximum. A positive 
correlation is found between the reconnection rate and the plasmoid velocity in case A 
(Figure [Hb). But in case B, there is no strong correlation between the reconnection rate and 
the plasmoid velocity (Figure EH). 

Figure [10] shows the magnetic field and velocity field in the case of no resistivity (case 
C, r/o = 0). The plasmoid begins to move upward slightly without reconnection, due to the 
non-equilibrium of the initial state (Figure [TOa). The rising plasmoid pulls up the magnetic 
field and stretches the current sheet, but reconnection does not start due to the lack of the 
resistivity. Therefore the plasmoid is still trapped by the closed magnetic field and stops 
rising (Figure [TUb). This result shows the importance of reconnection in plasmoid ejection. 
In fact, the reconnection based on numerical resistivity occurred in a later time, so the 
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plasmoid began to move upward again (Figure [TUB). 



4. Discussion 

Now we shall discuss the comparison of the simulation results with obse rvation. From 



soft X -ray observations of plas moid ejections in compact impulsive flares, IShibata et al. 



( 119951 ) and lShimizu et al. bund that there is a positi ve correlation betwe e n the plas- 



moid velocity and the rise velocity of the post-flare loop. iQiu fc Yurchyshynl (120051 ) also 
found a similar relation between the speed of coronal mass ejections (CMEs) and the recon- 
nection rate. The rise velocity of the post-flare loop is proportional to the reconnection inflow 
speed, so this means that ther e is a positive correlat i on be tween the plasmoid velocity and 



the reconnection inflow speed. iNarukage fc Shibatal (120061 1 discussed a correlation between 



CMEs speed and the inflow speed. This is actually obtained by the present result. 

How is this result explained physically? In case A, when the resistivity increases, the 
reconnection rate (t]J) increases. Since t]J = Vi n B in steady state, this means that the inflow 
speed (v in ) increases. Accordingly, the reconnection jet speed {v out ) also increases when v in 
increases because the mass conservations hold, v in L in = v out W out , in the incompressible limit. 

We also presented a physical explanation of how the plasmoid velocity is related to the 
enhancement of the resistivity in case B. In this case, we changed the plasmoid velocity by 
applying the additional force. We found weak dependence of the reconnection rate on the 
plasmoid velocity, on the contrary to the fact that the reconnection rate is supposed to be 
proportion to the plasmoid velocity because of mass conservation law. When the plasmoid 
velocity is larger, the plasmoid is accelerated farther from the X-point. When the plasmoid 
is situated far from the X-point, the plasmoid motion induces inflows only into the slow 
shock region near the plasmoid, but not into the X-point (Figure [TT1) . The dependence of 
reconnection rate on the plasmoid velocity is small. However the length of the slow shock 
increases when the plasmoid velocity is larger, so that energy release rate increases (Figure 
Eb). 

On the other hand there is a possibility that the positive correlation between the plas- 
moid velocity and the reconnection inflow speed holds in the initial phase, during which 
inflows will directly come into X-point. To examine the initial phase is however difficult in 
our model, because the initial configuration of our model is not in equilibrium and we cannot 
measure small changes of the reconnection rate and inflow speed in the initial phase. We 
will try to perform simulations with the initial condition in equilibrium in future. 



It should be mentioned that in our model (Figure [2]), magnetic flux stored in the plas- 
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moid is lost by the magnetic reconnection between the magnetic field overlying the plasmoid 
and magnetic flux stored in the plasmoid HI Therefore, if the amount of the magnetic flux 
stored in the plasmoid is small, the plasmoid disappears and does not propagate to the 
higher corona. There is another role of the reconnection occurring above an ejecting plas- 
moid. If the reconnection does not occur at the top of the plasmoid, the plasmoid might 
be decelerated by magnetic tension force, while if the reconnection occurs, the plasmoid is 
not decelerated. However in this simulation the plasma beta is high above the plasmoid (see 
Figure [2]), so the plasmoid is decelerated by gas pressure even if the reconnection occurs. 



5. Conclusion 



In this paper, we performed several MHD simulations to examine the basic physical 
relation between the plasmoid velocity and reconnection rate in the context of the plasmoid- 
induced-reconnection model for solar flares. The initial magnetic configuration, resistivity, 
and the additional force are parameters in these simulations. When we changed the am- 
plitude of resistivity (case A), the reconnection rate and the plasmoid velocity changed, 
showing a positive correlation. We showed that the reconnection rate (i.e. inflow speed) 
and the plas moid velocity are closely related to each other. This result is con sistent with 
observations (IShibata et al.lll995t iQiu fc Yurchyshynl 120051 ; IShimizu et al.ll2008l ) supporting 
the plasmoid-induced-reconnection model of impulsive flares. 



Numerical computations were carried out in part on general common use computer 
system at the Center for Computational Astrophysics, CfCA, of the National Astronomical 
Observatory of Japan. This work was supported by the Grant-in-Aid for the global COE 
program "The Next Generation of Physics, Spun from Universality and Emergence" from the 
Ministry of Education, Culture, Sports, Science and Technology (MEXT) of Japan and by 
the Grant-in-Aid for Creative Scientific Research of the MEXT "The Basic Study of Space 
Weather Prediction" (17GS0208, PI: K. Shibata). 



A. An acceleration mechanism of a plasmoid 

In this appendix, we examine an acceleration mechanism of a plasmoid semi-analytically. 
Let us consider the situation shown in Figure [12j In this situation, a plasmoid can be 
accelerated by the momentum added by the reconnection jet and the magnetic pressure 
gradient force, and decelerated by the magnetic tension force. Here, the magnetic tension 
force can be neglected in the magnetic field configuration in Figure [121 since there is almost 
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no overlying closed field above the plasmoid. (Note that the gravitational force is neglected 
since there is no gravity in our simulations, and the gas pressure can also be neglected in the 
region around the initial plasmoid.) In this magnetic field configuration, the force provided 
by the reconnection jet is estimated as 

Fjet = pvj et w ~ pv\w [dyn cm' 1 ], (Al) 

where p is the mass density in the jet, Vj e t is the velocity of the jet, which is comparable 
with the local Alfven velocity (va) in the inflow region, and w is the width of reconnection 
jet. The force exerted over the plasmoid due to is estimated as 

Fma 9 = W 2 ±f- ~ W 2 ±-f- [dyn cm" 1 ], (A2) 
ay 87r H 87T 

where W is the spatial size of the plasmoid and H(> W) is the spatial scale of the magnetic 
field. The ratio between these forces is as follows: 

F jet n (H\fw 



M = K^ = 2 \w) U < A3 > 

When M is greater than unity, the momentum provided by the reconnection jet is dominant. 
In this work, numerical simulation shows w/W ~ 1/4, H/W ~ 3/2 to 2 and M ~ 1 so that 
the these forces are comparable. 



REFERENCES 

Aschwanden, M. J. 2002, Space Science Reviews, 101, 1 
Carmichael, H. 1964, in The Physics of Solar Flares, 451 — h 
Chen, P. F., & Shibata, K. 2000, ApJ, 545, 524 
Coppi, B., & Friedland, A. B. 1971, ApJ, 169, 379 
Forbes, T. G. 1990, J. Geophys. Res., 95, 11919 
-. 2000, J. Geophys. Res., 105, 23153 

Forbes, T. G., & Malherbe, J. M. 1991, Sol. Phys., 135, 361 
Forbes, T. G., & Priest, E. R. 1983, Sol. Phys., 84, 169 
Hirayama, T. 1974, Sol. Phys., 34, 323 



-17- 



Hu, Y. Q. 1989, Journal of Computational Physics, 84, 441 

Isobe, H., Shibata, K., & Machida, S. 2002, Geophys. Res. Lett., 29, 10 

Karlicky, M. 2004, A&A, 417, 325 

Kim, Y.-H., Moon, Y.-J., Cho, K.-S., Bong, S.-C, & Park, Y. D. 2005, ApJ, 635, 1291 

Kliem, B., Karlicky, M., & Benz, A. O. 2000, A&A, 360, 715 

Klimchuk, J. 2001, in Space Weather, AGU Monograph 125, 143-157 

Kopp, R. A., & Pneuman, G. W. 1976, Sol. Phys., 50, 85 

Lin, J., & Forbes, T. G. 2000, J. Geophys. Res., 105, 2375 

Lin, J., Forbes, T. G., Isenberg, P. A., & Demoulin, P. 1998, ApJ, 504, 1006 

Lin, J., Soon, W., & Baliunas, S. L. 2003, New Astronomy Review, 47, 53 

Lin, J., & van Ballegooijen, A. A. 2005, ApJ, 629, 582 

Magara, T., Mineshige, S., Yokoyama, T., & Shibata, K. 1996, ApJ, 466, 1054 
Magara, T., Shibata, K., & Yokoyama, T. 1997, ApJ, 487, 437 

Masuda, S., Kosugi, T., Hara, H., Sakao, T., Shibata, K., & Tsuneta, S. 1995, PASJ, 47, 
677 

Masuda, S., Kosugi, T., Hara, H., Tsuneta, S., & Ogawara, Y. 1994, Nature, 371, 495 

McAllister, A., Uchida, Y., Tsuneta, S., Strong, K. T., Acton, L. W., Hiei, E., Bruner, M. E., 
Watanabe, T., & Shibata, K. 1992, PASJ, 44, L205 

Narukage, N., & Shibata, K. 2006, ApJ, 637, 1122 

Nitta, N., & Akiyama, S. 1999, ApJ, 525, L57 

Ogawara, Y., Takano, T., Kato, T., Kosugi, T., Tsuneta, S., Watanabe, T., Kondo, I., & 
Uchida, Y. 1991, Sol. Phys., 136, 1 

Ohyama, M., & Shibata, K. 1997, PASJ, 49, 249 

-. 1998, ApJ, 499, 934 

Parker, E. N. 1963, ApJS, 8, 177 



-18- 



Petschek, H. E. 1964, in The Physics of Solar Flares, 425-+ 
Priest, E. R., & Forbes, T. G. 1986, J. Geophys. Res., 91, 5579 
-. 2002, A&A Rev., 10, 313 

Qiu, J., & Yurchyshyn, V. B. 2005, ApJ, 634, L121 

Sakajiri, T., Brooks, D. H., Yamamoto, T., Shiota, D., Isobe, H., Akiyama, S., Ueno, S., 
Kitai, R., & Shibata, K. 2004, ApJ, 616, 578 

Shafranov, V. D. 1966, Reviews of Plasma Physics, 2, 103 

Shibata, K. 1996, Advances in Space Research, 17, 9 

Shibata, K. 1997, in ESA SP-404: Fifth SOHO Workshop: The Corona and Solar Wind 
Near Minimum Activity, 103 — h 

-. 1999, Ap&SS, 264, 129 

Shibata, K., Ishido, Y., Acton, L. W., Strong, K. T., Hirayama, T., Uchida, Y., McAllister, 
A. H., Matsumoto, R., Tsuneta, S., Shimizu, T., Hara, H., Sakurai, T., Ichimoto, K., 
Nishino, Y., & Ogawara, Y. 1992, PASJ, 44, L173 

Shibata, K., Masuda, S., Shimojo, M., Hara, H., Yokoyama, T., Tsuneta, S., Kosugi, T., & 
Ogawara, Y. 1995, ApJ, 451, L83+ 

Shibata, K., & Tanuma, S. 2001, Earth, Planets, and Space, 53, 473 

Shimizu, M., Nishida, K., Takasaki, H., Shiota, D., Magara, T., & Shibata, K. 2008, ApJ, 
683, L203 

Shimojo, M., Hashimoto, S., Shibata, K., Hirayama, T., Hudson, H. S., & Acton, L. W. 
1996, PASJ, 48, 123 

Shiota, D., Isobe, H., Chen, P. F., Yamamoto, T. T., Sakajiri, T., & Shibata, K. 2005, ApJ, 
634, 663 

Sterling, A. C, & Moore, R. L. 2005, ApJ, 630, 1148 

Sturrock, P. A. 1966, Nature, 211, 697 

Sui, L., & Holman, G. D. 2003, ApJ, 596, L251 



-19- 

Sweet, P. A. 1958, in IAU Symposium, Vol. 6, Electromagnetic Phenomena in Cosmical 
Physics, ed. B. Lehnert, 123 — h 

Takasaki, H., Asai, A., Kiyohara, J., Shimojo, M., Terasawa, T., Takei, Y., & Shibata, K. 
2004, ApJ, 613, 592 

Tsuneta, S., Hara, H., Shimizu, T., Acton, L. W., Strong, K. T., Hudson, H. S., & Ogawara, 
Y. 1992a, PASJ, 44, L63 

Tsuneta, S., Takahashi, T., Acton, L. W., Bruner, M. E., Harvey, K. L., & Ogawara, Y. 
1992b, PASJ, 44, L211 

Ugai, M. 1996, Physics of Plasmas, 3, 4172 

Ugai, M., & Tsuda, T. 1977, Journal of Plasma Physics, 17, 337 

Yokoyama, T., & Shibata, K. 2001, ApJ, 549, 1160 

Zhang, J., Dere, K. P., Howard, R. A., Kundu, M. R., & White, S. M. 2001, ApJ, 559, 452 



This preprint was prepared with the AAS IATgX macros v5.2. 



-20- 



(a)t=40 (b)t=45 (c)t=50 (d) t=55 




Fig. 4. — The magnetic field and velocity field nearby the plasmoid in case A, 770 = 0.002. 
The solid lines denote magnetic field lines, and arrows represent velocity field. Magnetic 
reconnection occurs at the top of the plasmoid. 
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(a) Drift velocity 
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(c) Plasmoid velocity 
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(b) Reconnection rate 
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Fig. 5. — These panels show the temporal variation of the standard case (case A, r] = 0.002). 
(a) Maximum drift velocity of = j z /p in the current sheet. Dotted line shows — v c . (b) 
Maximum reconnection rate in the current sheet, not normalized (i.e. dip/dt = rjj). (c) 
Plasmoid velocity, (d) Inflow speed measured at the point of maximum reconnection rate. 




Fig. 6. — Distribution of ^-component of current density (j z ) at t = 30 in the standard case 
(case A, rj = 0.002). The solid lines denote magnetic field lines. 
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Fig. 7. — Panel (a) shows the initial configuration of magnetic field, in which solid lines 
represent the distribution of magnetic flux ip. At panels (b) and (c) the contours show the 
time evolution of the point on each field line that is located at the y-axis in the standard 
case (case A, 770 = 0.002). The colors variation show the temporal change of gas pressure 
in panel (b) and vertical velocity v y in panel (c). The dotted lines in panel (b) show the 
trajectory of the fast shocks (termination shocks) formed at the top of the flare loop and 
the bottom of the plasmoid. The white line in panel (c) shows the position of the upward 
jet where v y > v . 
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(c) Plasmoid velocity vs. Reconnection rate 
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Fig. 8. — Time variation of (a) the reconnection rate (dip/dt) and (b) the plasmoid velocity 
(vpiasmoid) in case A, when free parameter r] is changed. Panel (c) shows the correlation 
between the maximum value of the reconnection rate and the plasmoid velocity when the 
reconnection rate becomes the maximum in case A. The black diamond is the standard case, 
7/o = 0.0002. 
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(a) Reconnection rate 
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(c) Plasmoid velocity 
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(b) Energy release rate 



Upward +1.0 
Upward +0.5 
No force 
Downward -0.5 
Downward -1.0 




(d) Plasmoid velocity vs. Reconnection rate 
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Fig. 9. — Time variation of (a) the reconnection rate (dip/dt), (b) the energy release rate, 
and (c) the plasmoid velocity (v p i asrnoi d) in case B, when the additional force F y is changed. 
Panel (d) shows the correlation between the maximum value of the reconnection rate and the 
plasmoid velocity when the reconnection rate becomes the maximum. The white diamonds 
show the result in in case B. Solid lines in panel (a)(b)(c) and the black diamond in panel 
(d) show the result in standard case (case A, i]q = 0.002), just for comparison. 
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(a)t=15 (b)t=30 (c)t=60 (d) t=90 




Fig. 10. — The magnetic field and velocity field in the case of no resistivity (case C). The 
solid lines denote magnetic field lines, and arrows represent velocity field. 



(a) Upward +1 



(b) No Force (c) Downward -1 




-0.2 0.0 0.2 0.4 0.6 0.8 

X 



Velocity: 



Fig. 11. — The magnetic field and velocity field nearby the inflow region in case B, when 
the additional force F y is changed. The solid lines denote magnetic field lines, and arrows 
represent velocity field. The times are t = 14, 18, 28 respectively, when the position of the 
plasmoid is y — 1.5 approximately. 
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Fig. 12. — Schematic illustration of an acceleration mechanism of a plasmoid. The magnetic 
field configuration is the result of the simulation at t = 20. H(> W) is the spatial scale of 
the magnetic field, which can be consider as the distance between the inflow region and the 
top of the plasmoid where the magnetic pressure is low. 



